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FOREWORD 


This  report  presents  a  simple  prescription  for  generating  the  appropriate  contour  mtegr^s 
describing  the  Green’s  function  for  propagation  in  a  waveguide.  It  then  proceeds  to  descnbe 
methods  by  which  these  contour  integrals  can  be  evaluated  in  terms  of  sums  over  normal  modes 
and  branch  cut  integrals.  Examples  for  the  one-  and  two-layered  waveguide  are  developed  in 

detail. 
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INTRODUCTION 

This  report  presents  a  simple  prescription  for  generating  the  appropriate  contour  integrals 
describing  the  Green's  function  for  propagation  in  a  waveguide.  It  also  describes  a  method  by 
which  these  contour  integrals  can  be  evaluated  in  terms  of  sums  over  normal  modes  and  branch 
cut  integrals. 


INTEGRAL  REPRESENTATIONS  AND  MULTIPATHS 

The  first  assumption  is  that  the  source  is  a  harmonic  point  source  with  a  time  dependence 
of  the  form  .  Second,  the  Green  function  for  a  point  source  in  free  space  is  given  by  the 
expression 

G(r, /)  =  (1) 

which  can  be  represented  in  the  form 

G(r,  =  i/SnjZ  (2) 

in  terms  of  the  cylindrical  functions  appropriate  for  a  waveguide.  Here,  q  is  the  horizontal 
wavenumber  and  h  is  the  vertical  wavenumber  of  the  field.  The  conventions  of  Ewing,  Jardetsky 
and  Press  for  the  branch  cut  of 

=  (3) 

is  chosen  such  that  the  imaginary  part  of  h  is  positive  on  the  physical  sheet.  For  real  frequencies, 
the  branch  cut  consists  of  a  cut  from  +k  to  the  origin  to  +ioo.  The  cut  for  negative  frequencies  is 
obtained  by  the  making  the  transformation  h-^—h  along  the  cut.  The  z-axis  is  oriented  towards 
the  bottom  of  the  waveguide.  This  convention  affects  the  signs  in  the  exponential  of  the 
cylindrical  form  of  the  Green  function. 

Using  Eqn  3  for  the  Green  function  one  can  obtain  an  integral  representation  of  the  Green 
function  in  a  waveguide.  Consider  a  source  and  receiver  in  a  homogenous  layer,  where  U  and  V 
are  the  total  reflection  coefficients  of  the  upper  and  lower  boundaries  of  this  layer.  Then  the 
Green  function  for  the  wave  guide  may  be  expressed  in  the  form 

G^ry)  =  i/Sirn  V9-P)-F(z,z0  (4) 

where  the  function  F(z,z')  is  given  by  the  following  multiple  scattering  summation  of  multipaths. 
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^  +  U_e-iK2a-z-z')  +  y  gih(2b-z-z')  +  (/V.e2‘W-</..(z-z')).((7V.e^‘'“')"  (5) 

This  summation  consists  of  the  contributions  due  to  the  direct  path,  a  surface  bounce,  a  bottom 
bounce,  and  a  bottom  bounce  followed  by  a  surface  bounce  and  all  multiple  bounces.  Here,  a  is 
the  z  coordinate  of  the  upper  surface  and  b  is  the  z  coordinate  of  the  lower  surface  and  d  =  b  -  a  is 
the  depth  of  the  layer.  In  deriving  the  above  expression  it  was  assumed  that  z  >  z'.  The 
expression  for  z  <  z'  is  obtained  by  switching  z  and  z'  in  the  above  expression. 

Eqn  5  can  be  further  simplified  in  the  form 

F(z,  z')  =  1  -  UVe^^'^)  (6) 

appropriate  for  a  contour  integral,  where  z  >  z'.  The  case  of  z  <  z'  is  represented  by  switching  z 
and  z'  in  the  above  expression. 

The  contour  integral  representation  of  the  waveguide's  Green  function  is  given  by  the 
following  expression. 

G.(r,  F)  =  i/Sn.  (7) 

Consider  the  generalization  where  the  source  and  the  receiver  are  in  different  layers.  Let 
the  source  z'  be  in  layer  1,  where  a'  and  b’  are  the  z  coordinates  of  the  upper  and  lower  surfaces. 
Let  the  receiver  z  be  in  layer  2,  where  a  and  b  are  the  z  coordinates  of  the  upper  and  lower 
surfaces  of  this  layer.  Then  using  the  multipath  expansion  of  the  waveguide  one  can  express  the 
Green  function  for  the  waveguide  in  the  following  form, 

G,(r,r')  =  im4Z  *.Hi‘’(«P).F.(z').H'l2.f2(z)  (8) 

where  W12  is  the  total  transmission  coefficient  between  layer  1  to  layer  2.  Assuming 
z  >  z',  the  function  Fi  (zO  is  obtained  by  summing  all  multipaths  which  strike  the  interface  b'.  The 
function  F2(z)  is  obtained  by  summing  all  multipaths  from  the  upper  surface  a  to  the  receiver  at  z. 
In  the  case  z  >  z',  these  functions  are  given  by  the  following  expressions. 

Fi(zO  =  +  ’  (9) 

F2(z)  =  - 1/2^2. (10) 

In  the  above  expression  Ui  is  the  total  reflection  coefficients  at  a’  of  all  the  layers  above  a',  vi  is 
the  surface  reflection  coefficient  at  b'  as  if  the  layer  below  b'  were  a  homogeneous  halfspace,  U2  is 
the  total  reflection  coefficient  at  a  of  all  the  layers  above  a,  and  V2  is  the  total  reflection 
coefficient  at  b  of  all  the  layers  below  b.  Note  the  use  of  the  vertical  wavenumber  in  layer  1  in  the 
function  Fi  (zO .  and  the  use  of  the  vertical  wavenumber  in  layer  2  in  the  function  F2(z) . 
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In  the  case  z  <  z\  the  function  F\  {z')  is  the  sum  of  all  multipaths  from  the  source  to  the 
upper  interface  at  a',  and  the  function  F2(z)  is  the  sum  of  all  multipaths  from  the  lower  interface  at 
b  to  the  receiver.  In  this  case,  these  functions  are  given  by  the  following  expressions, 

F, iz')  =  +  Vi  -^'))/(l  -  Ml .  Vi  )  (11) 

F2(z)  =  (12) 

where  u\  is  the  surface  reflection  coefficient  at  a'  as  if  the  layer  above  a'  were  a  homogeneous 
halfspace,  and  Vi  is  the  total  reflection  coefficient  at  b'  of  all  the  layers  below  b'. 

Note,  this  form  of  the  integral  representation  is  only  valid  in  the  case  both  the  source  and 
the  field  point  are  in  fluid  layers.  However,  the  intervening  layers  may  be  elastic  or  Biot  models, 
since  this  form  of  the  integral  representation  needs  only  the  reflection  and  transmission  coefficient 
of  these  intervening  layers.  An  integral  representation  for  an  elastic  or  Biot  waveguide  would 
require  use  of  the  propagator  matrix  formalism  to  represent  the  field. 

NORMAL  MODES 

This  section  describes  the  normal  mode  representation  one  can  obtain  from  the  previous 
integral  representations  of  the  Green  function  by  representing  these  contour  integrals  in  terms  of  a 
sum  of  residues  from  poles  and  branch  cut  integrals.  For  the  sake  of  simplicity,  only  the  case 
where  both  the  source  and  the  receiver  are  in  the  same  layer  is  treated. 

The  Green  function  for  a  source  and  receiver  in  the  same  layer  is  given  by  the  following 
contour  integral  representations  in  the  complex  h-plane. 

Gw(r,  F)  =  t'/STt.  +  t/.e^‘^<).(e^'^>  +  V.e2'^.e-'^>)/(l  -  (13) 

This  integral  may  be  represented  by  the  contour  integral  in  the  q-plane  as 

G^(r,r^)  =  i787t.Jl^if^//i./fo  ^(<?p).(e-'^  +  (14) 

In  general  the  poles  of  the  integrand  lie  on  or  near  the  branch  cuts.  In  general,  this  contour 
integral  may  be  represented  as  a  sum  of  pole  contributions  and  a  cut  contributions.  The  cut 
contributions  occur  when  the  integrand  is  not  an  even  function  of  h.  In  the  case  U  =  -1  and  V  = 
+1  for  a  pressure  release  surface  and  a  rigid  bottom,  the  integrand  is  an  even  function  of  h,  hence 
the  cut  contribution  is  zero.  In  the  case  U  and  V  have  the  symmetry  U(-h)  =  U(h)  '  and  V(-h)  = 
V(h)''  ,  the  integrand  is  an  even  function  of  h  and  its  branch  cut  contribution  vanishes.  In  general, 
the  cut  contribution  due  to  any  bounded  layer  is  zero.  Cut  contributions  generally  occur  when 
you  have  an  unbounded  region  at  the  ends  of  the  waveguide. 

The  branch  cut  for  q  =  in  the  complex  h-plane  is  chosen  such  that  Im(q)  >  0  in 

the  complex  h-plane.  The  above  contour  can  be  completed  by  adding  the  contour  at  infinity  in  the 
upper  half  plane,  since  the  integrand  is  exponential  damped  along  the  added  contour.  In  addition 
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to  the  branch  cut  for  q,  there  exists  branch  cuts  for  the  vertical  wavenumber  in  the  other  layers 
and  if  there  exists  an  unbounded  region  in  the  waveguide  there  will  be  a  cut  contribution  to  the’ 
above  integrand  in  addition  to  the  pole  contributions. 

The  poles  of  the  above  integrand  correspond  to  the  complex  solutions  of  the  equation. 


(1  -  =  0 


(15) 


The  pole  contributions  to  the  above  integrand  are  of  the  form 

where 


A'  =  |(1  -  =  -2i.{d+ 1  ^.|(t7V)) 


(17) 


The  normal  mode  representation  of  the  Green  function  is  simply  the  sum  over  all  the 
residues  from  the  poles  and  cut  contributions  in  the  upper  half  q-plane. 

HOMOGENEOUS  WAVEGUIDE  WITH  A  RIGID  BOTTOM 

This  section  derives  the  Green  s  function  for  a  homogeneous  waveguide  with  a  pressure 
release  surface  and  a  rigid  bottom.  The  reflection  coefficients  U  and  V  are  given  by  the 
expressions  U  =  -1  and  V  =  +1.  Substituting  these  values  into  Eq.  16  for  the  pole  contributions  to 
the  Green  function  the  following  expression  is  obtained, 

^(9P).(e'^'^'  -  -  g-^z) 

which  may  be  represented  in  the  form 

where 

sin(fe) 

are  the  normalized  depth  function  of  the  waveguide. 

The  characteristic  equation  for  the  waveguide  takes  the  form 
(l+e2'^)  =  0 

whose  solutions  are  given  by  the  expression 
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/z  =  (n-l/2)*7l/i  (21) 

forn=  1,2... 

TWO  LAYERED  WAVEGUIDE  WITH  RIGID  BOTTOM 

This  section  applies  the  previous  sections  in  the  case  the  waveguide  consists  of  2  fluid 
layers,  where  the  top  layer  has  a  pressure  release  surface  and  the  bottom  layer  has  a  rigid  bottom. 
In  this  case  the  reflection  coefficients  at  the  surface  and  bottom  are  given  by  the  expressions. 

=-l  (22) 

V2=V2=+1  (23) 

The  reflection  coefficients  at  the  interface  of  these  2  layers  is  more  complicated.  Note  that  the 
surface  reflection  coefficients  at  this  interface  are  given  by  the  expressions, 

v\  =  {m.h]  - h2)lim.h\  +  h2)  (24) 

W2=-V)  (25) 

where  m  is  the  ratio  of  the  density  in  layer  2  and  the  density  in  layer  1.  The  total  reflection 
coefficients  at  this  interface  are  given  by  the  expressions 

V,  =  (vi  +  +  vi  V2.e^^''‘2‘'')  =  (vi  +  e+2<A2^2)/(i  +  vi  (26) 

U2  =  (m2  +  Ui.e+2‘^‘^')/(l  +M2Ui.e+2‘^‘'i)  =  (-vi  +  (27) 

In  the  case  both  the  source  and  receiver  are  in  the  upper  layer  (layer  1),  the  Green  function 
is  given  by  the  following  contour  integral. 

G^{r,  /)  =  iVSTtH  77 +  Vi +  Vi (28) 
The  normal  mode  (pole)  contribution  to  this  integral  is  given  by  the  expression, 
-l/4.;^.H|)‘\^p).(e-^''"^'  -  e-‘'>i^')(e+'''i^  -  (29) 

where 

A'=-2/(rfi+^.|;(V,))  (30) 

and  the  normal  modes  (poles)  are  solutions  of  the  following  equation. 


5 


CSSArR-95/10 


A  =  (l  +  Vi 


(31) 


Multiplying  the  characteristic  equation  by  {mh\  +/z2)(l  +  one  can 

express  the  poles  (normal  modes)  as  solutions  of  the  equation. 

T  —  2.(mh]  —  h2)cos(h\d\  —  h2d2)  +  2.{mh\  +/Z2)  *  cos(h\d]  +/i2<5?2)  (32) 

Note,  this  function  is  a  real  odd  function  of  the  vertical  wavenumber  in  both  the  first  and  second 
layer,  which  implies  that  if  Aj  =  a  is  a  solution  of  this  equation,  then— a,  a*,— a*  are  also 
solutions  of  the  characteristic  equation.  The  roots  to  this  equation  are  real  for  real  wavenumbers 
k\,k2.  Hence  the  roots  lie  along  the  branch  cuts  for  the  vertical  wavenumber. 


In  order  to  compute  the  derivative  of  the  characteristic  equation  define  the  quantity 


f=k]-ki 


(33) 


The  vertical  wavenumber  in  the  second  layer  maybe  rewritten  in  the  form 

hi  =  (■ 

in  terms  of  the  wavenumber  in  the  first  layer.  The  derivative  of  the  characteristic  equation 
evaluated  at  a  solution  of  the  characteristic  equation  is  given  by  the  expression. 

=  -2idi  -■^^{2id2.h\,i{mhi  +h2)  +  {mhi  -h2).e'^'^‘'''‘‘')  +  ({mh2+hi)  +  {mh2-hi).e*^‘''^‘‘^) 
Himh2-hi)  +  imh2+hi).e*^‘^^‘‘^).e*^‘>'>d^}/((mhi  +/i2)  +  (/n/ii -h2).e*^‘''^‘‘2)  (3 

In  the  homogeneous  limit,  m  =  1 ,  h2=h]  =  h,the  characteristic  equation  and  its  derivative  are 
given  by  the  expressions 


A  =  (1  +  =  0 


(36) 


and 

A^  =  — 2i.(rfi  +£^2) 

appropriate  for  a  homogeneous  waveguide  with  a  rigid  bottom. 


(37) 
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In  the  case  the  source  z'  is  in  the  top  layer,  and  the  field  point  z  is  in  the  bottom  layer,  the 
Green  function  is  given  by  the  following  expression, 

G..(r,  /)  =  mn.  IZ  iz').Wn.F2iz)  (38) 

where 

W^i2  =  (l+v,)  (39) 

is  the  transmission  coefficient  between  layers  1  and  2,  and 

F|  (zO  =  (40) 

and 

F2{z)  =  (e'^‘^U-:2)  +  g+2-V.2<i2_g-rt2U-Z2))/(l  _  (41) 

Here,  z\  =  0,Z2  =  +<^i,Z3  =  {d\  +d2)  are  the  z  coordinates  of  the  three  interfaces  of  the 
waveguide.  By  making  use  of  the  identity, 

(1  +vi.e^2'''>‘^‘).(l  -G2.e+2‘'''2‘^2)  =  (1  +v,.e+2<'>2rf2).(n.  v,.e+2<V!,rfi)  (42) 

the  expression  for  the  Green  function  can  be  expressed 
in  the  following  form. 

G„.(r, F)  =  -l/lK. IZ  +  V, )/(l  +  vi .e+2.7.2<i2)}.e+-v.,</,+iM2/A 

sin(/iiz').cos(/i2(z-rfi -c?2))  (43) 

The  normal  modes  are  solutions  of  the  characteristic  equation 

A  =  (l  +  Vi.e+2//.,</,)^0  (44) 

and  the  pole  contribution  to  this  integral  is  of  the  form. 

(+i/A').i{il+V])/i\+vi.e'^^‘^^‘^^)}.e'^'''‘‘''*'‘^^‘‘^.HQ\qp).smih\z').cosih2{z-d\  -d2))  (45) 

This  solution  can  be  recast  into  the  form 

ilA.H^^\qm^{z')Mz)  (46) 

where  f\{z)  and  /2(z)  are  the  functions 
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f\(z)  =  J-4i/A' .  sin(/iiz)  (47) 

and 

/2(z)  =  J^.(+0{(l  +  vi)/(l  g+iAi</i+tft2<i2,cos(/z2.(z-<3fi -dj)).  (48) 

Similarly,  the  solution  in  the  case  both  the  source  and  field  point  are  contained  within  layer  1  can 
be  recast  into  the  form. 

i/4.Ho\^p)./i(zO/iW  (49) 

In  the  homogeneous  limit,  m  =  1,  h\  =  h2-h,i\\&  depth  functions  are  equal  to  the 
following  expression, 

A  (z)  =/2(z)  =  .  sm(hz)  (50) 

appropriate  for  a  homogeneous  waveguide. 

FLUID  LAYER  OVER  A  HOMOGENEOUS  HALFSPACE 

This  section  will  describe  the  Green's  function  for  a  fluid  layer  over  a  fluid  halfspace.  In 
this  case,  the  reflection  coefficients  are 

Ui=M,  =-l  (51) 

Vi  =V]  =(mhj -h2)/(mh] +h2)  (52) 

U2=-Vi  (53) 

f/2  =-(v]  ■/>,</,)  (54) 

V2  =  V2  =  0  (55) 

where  the  conventions  of  the  previous  section  have  been  adopted. 

In  the  case  z  >  z',  and  the  source  and  the  receiver  are  contained  within  layer  1,  the  Green's 
function  is  given  by  the  contour  integral. 

GUry)  =  i/87t. +  +  (56) 
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The  normal  modes  are  solutions  of  the  characteristic  equation, 

A  =  (l =  0 

which  can  be  recast  into  the  form. 

A  =  ((m/ll  +  hi)  +  (m/ll  +/12)  =  C  (58) 

The  normal  mode  contribution  to  this  contour  integral  is  given  by  the  residue 
^.//o\^p).sin(/iizO-sin(/iiz),  (59) 

where 

A'  =-2id]  -•^gly(vi)  =-2ii/i  -  —.{(mhi  - h\)l{mh  \  +  hi)  -  {mhi  +  h\)l{mh  \  +/12)}  (60) 

is  the  derivative  of  the  characteristic  equation  evaluated  at  the  pole. 

The  integrand  is  an  even  function  of  Ai ,  but  it  is  not  an  even  function  of  hi 
This  leads  to  a  discontinuity  across  the  hi  cut  in  the  complex  q-plane,  which  leads  to  a  branch  cut 
contribution  to  the  above  contour  integral,  in  addition  to  the  pole  contributions.  This  branch  cut 
contribution  is  due  to  the  evanescent  fields  generated  by  the  halfspace,  and  it  is  of  the  form, 

j;~  hi.dhi.(F{+h, ,  +/12)  -  F(-/ii,  -  hi))  (61) 

along  the  branch  cut  of  hi  and  F{h\hi)  is  the  function. 

F{hx,  hi)  =  (i78Tt)f^//o  \^p).(c-''"^'  -  e+''’'^')(e+''"^  +  vi  .e-*iz)/(l  +  vi  (62) 

Making  use  of  the  identities, 

Vi( — /12)  — +l/vi  (+/12)  (63) 

v^-Zii)  =+l/vi(+/ii)  .  (64) 

this  contour  integral  can  be  written  in  the  form. 

-i/27C.  d/i2./z2//ii .^0  \9P).  sin(/iiz0- sin(/iiz).(v]  -  l)/((vi  + 

=  +i/271.  j^d/i2.Ho\^p).sin(/iizO-sin(/iiz).m./i|/((m/iicos(/ii(ii))^  +  (/i2.sin(/iidi))^)  (65) 
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The  integrand  is  an  even  function  of  h\  and  h2 . 


Assuming  k2^»  I  and  using  the  relationship 

q  =  +ijhl-kl  (66) 

the  integrand  is  exponentially  damped  for  h2  greater  than  k2  due  to  the  exponential  behavior  of 
the  Hankel  function  for  large  arguments.  In  this  case,  the  chief  contribution  to  the  integrand 
comes  from  the  oscillatory  portion  of  the  integrand  from  h2  =  0  to  k2. 

Assuming  /:2P  »  1  i  the  branch  line  integral  can  be  re-expressed  in  the  form, 

jj”  hldh2.e^‘^P.F(h2)  j.r(3/2).(2k2/p)^^^.e*‘^^P-‘-^”''^.F(h2  =  0)  (67) 

where  the  asymptotic  expansion  ( Eq.  68  )  has  been  used 

^o^(^P)->  (68) 

for  the  Hankel  functions,  and  the  method  of  steepest  descent  to  evaluate  integrals  of  the  form 

j^x'’.dx.e^'^^\F(x)  r((n+  l)/2).l/'(xo)/2|-(”+')^.e+'^'>>^(''^')'^^F(xo)  (69) 

where  xo  is  the  saddle  point  satisfying  the  equation,  f(xo)  =  0.  The  asymptotic  expression  for  the 
branch  line  integral  is  given  by  the  expression, 

-i.k2/(2nm).e'*'‘^^<^/(yp)^ .  sin(yz').  sin(yz)/cos('Ydi)^  (70) 

where  Y=  Jk^  -k\  .  Note  in  particular,  that  the  asymptotic  form  of  the  branch  line  integral  falls 
off  as  p~^ ,  whereas  the  asymptotic  form  of  the  normal  mode  contributions  fall  off  as  p"'^ .  Thus, 
the  branch  line  integral  is  negligible  at  far  ranges,  where  the  normal  mode  terms  dominate. 

Let  us  next  consider  the  case  the  source  is  in  layer  1  and  the  field  point  is  in  layer  2.  In 
this  case,  the  Green's  function  is  given  by  the  expression 

The  normal  mode  contribution  to  this  integral  is  of  the  form 
-i/27i.(H-V]).e+'''“'‘/A'./7[,'\<?p).sin(/n(z')-e^^^^^“‘''L  (72) 
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The  branch  cut  integral  is  of  the  form 

Uln.  h2dh2lh\.HQ\qp).s\n{h\z'){\+v\).{v\.sm{h2{z-d\)-h\d\)-¥sm{h2{z-d\)  +  h\d\)) 

g+2l/ll<<l)  (yj  (73) 

This  can  be  rewritten  in  the  form 

UlTl.  h2.dh2.HQ\qii).sm{h\z').{mhi  .sm{h2{z  -  d\)).cos{h\d\)  +  h2COs{h2{z-  dx))s,\n{h\d\)} 
.m/{(m/iicos(/iitfi))^  +  (/i2sin(/iiffi))^)}.  (74) 

This  integral  is  of  the  form 

hl.dh2.e^'^^ .FQ12)  (75) 

where 

F(/i2)  =  Uln.  sm(h\z').{mh]/h2.  sm{h2iz-d\)).cosih\d\)+ 

cos{h2.(z  - d])).  sm{h]dx)}/{{mh]COs{h\d\))^  +  ih2sm{h\d\))^)}  (76) 

An  asymptotic  expression  for  this  branch  cut  integral  may  be  obtained  as  in  the  previous  case. 

The  resulting  asymptotic  expression  is  of  the  form 

-i./:2/(2ji.m).e'^'^P/(yp)2.  sin(YzO{sin(Trfi)  +  myXz-di).  cos(yJi)}/  cos(Y^fl  )^ .  (77) 

Note,  this  term  changes  linearly  with  respect  to  the  depth  of  the  field  point. 
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